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Abstract 

We derive the Hamilton equations of motion for a constrained system in the 



> 

form given by Dirac, by a limiting procedure, starting from the Lagrangean 
for an unconstrained system. We thereby ellucidate the role played by the 

O 

£T) ' primary constraints and their persistance in time. 

O ' 

o ■ The Hamiltonian formulation for systems whose dynamics is described 

by a Lagrangean with singular Hessian has been given a long time ago by 
Dirac [1], and has been elaborated ever since in numerous papers 2 . As is well 
^ ■ known all gauge theories fall into the class of singular systems. The usual 

starting point for deriving the Hamilton equations of motion is the singular 
Lagrangean which leads, in the language of Dirac, to primary constraints. 
While the primary constraints have no analog on the Lagrangean level, re- 
quiring their persistance in time leads to equations relating the coordinates 
and velocities. These are part of the Euler-Lagrange equations of motion. 
This persistance requirement may lead to secondary constraints and is the 
first step in the Dirac algorithm generating the constraints. The purpose of 
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2 It is impossible to quote here the very large number or papers. Comprehensive dis- 
cussions can be found e.g. in [2] 
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this paper is to derive the Hamilton equations of motion for systems with a 
singular Hessian, by starting from an unconstrained system and taking an 
appropriate limit. The formulation of the equations of motion within an 
extended phase space is quite natural, and the role played by the primary 
constraints, and the demand for their persistance in time will be illuminated 
thereby. 

Consider the following Lagrangean quadratic in the velocities 

L = \ H w ij(<l> a )^j - 2 Viiv)* - v (<l) • (!) 

ij i 

where q stands for the set of coordinates {<&} and a stands collectively for 
a set of parameters. For a ^ a c we assume that det W ^ 0, so that we are 
dealing with a non-singular system. We assume that the singular system of 
interest is realized for a = a c (the subscript c stands for "critical") where 
det W(q; a c ) = 0. Our aim is to derive the Hamilton equations of motion 
for the singular system, by taking the limit a — > a c of the Hamilton equa- 
tions of motion following from (l). 3 For a ^ a c , one readily constructs the 
Hamiltonian 

H = \ Efe + rfiW^fa + Vj ) + V{q) , (2) 
where the canonical momenta are related to the velocities by 

Pi = J2 w *Mi a )<lj ~ ■ (3) 

3 

The symmetric matrix W can be diagonalized by an orthogonal transfor- 
mation, Wd = C T WC with Ca = vf\ where are the orthonormalized 
eigenvectors of W , 

W(q; a)^\q; a) = X e (q; a)^(q; a) . (4) 

In terms of the eigenvalues and eigenvectors, W can be written in the form 

H/,-=EW\ (£) . (5) 
i 

3 Clearly the form of (1) is not unique. 
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Correspondingly W^ 1 is obtained by making the replacement A^ — > j^. Hence 
the Hamiltonian (2) is given by 

H=\E^ + V(q), (6) 

where 

Mq,p; «) := (p + ifa)) • «) • (7) 

For a ^ a c the Hamilton equations of motion then take the form 

OH ^ 1 <9<^ 
% ^ A £ % 

Consider now the limit a — > a c where cfet = 0. Let {A^ } denote the set 
of eigenvalues which vanish in this limit. In order to implement the limit we 
first write (8) and (9) in the form 

Pl = "2 j£ } % U^J - 2 § % U> J - ^7 • (ID 
Now, from (3) and (7) we have that 

^ = Y.vf 0) W l3 q 3 = \ k) ^-'q. (12) 

i 

Hence 

lim — £o = p €o (g, g; a c ) , (13) 

where 

Pi (q,q;a c )='q-v (eo \q;a c ) . (14) 
Note that the finiteness of the velocities in (12) implies that 

(f>to(q,P,at c ) = , (15) 
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which are just the primary constraints. Hence for a — > a c (10) reduces to 

f ^)V A » *</„, V V a Pi /„=„„ 

Consider next eqs. (11). Making again use of (12) and the fact that -r^ L \ a =a c - 
one obtains that for a — > a r these reduce to 



5 ^ ^2 



where use has been made of (13). 

We now notice that the first sum on the RHS of eqs. (16) and (17) are 
just given by 7^ and — ttS where H c is the canonical Hamiltonian obtained 
from (6) by taking the limit a — > a c : 

H c = \H yt^-t^^p; Oc) + (18) 

Here use has again be made of (15) and of the fact that lim Q ,^ Q , c ^ is finite. 
Hq is just the canonical Hamiltonian derived from the (1) for a = a c , eval- 
uated on the primary surface. Hence the equations of motion (16) and (17) 
take the well known form 

* = (w) 

dpi 

Pi = (20) 
where Ht is the "total" Hamiltonian 

H T = H c + Y,Pe 4>e , (21) 

to 

and /9^ are the undetermined projections of the velocities on the zero modes 
(14). Note that the derivatives in (19) and (20) are understood not to act 
on pt . Eqs. (19) and (20) must be supplemented by the primary constraints 
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(15) . In fact these equations only have a solution if qi(t) and Pi(t) are points 
in the submanifold defined by the primary constraints [3]. Note that from 

(16) it follows that the projection of q onto the zero mode reduces to 
an identity, since ^Zivf ' = 0. Hence we have obtained the Hamilton 
equations of motion for a constrained system in the form given by Dirac, 
by taking the limit a — > a c of the equations of motion for an unconstrained 
system. The primary constraints are just the statement that the projected 
velocities (14) are finite in this limit. 

Actually, eqs. (19) and (20), together with the primary constraints (15) 
do not directly yield the complete set of Lagrange equations of motion. 
These follow by also implementing the persistance in time of the primary 
constraints. The primary constraints themselves have no analog on the La- 
grangean level but allow us to recover the connection between momenta and 
velocities needed to express the Hamilton equations of motion in terms of 
Lagrangean variables. Thus from (16) and (7) it follows that 

ft= E yfa? + EOT' ^°Vf 0) > (22) 

where it is understood that we have set a = a c . Define the matrix 
constructed from and X e (£ ^ {£ }): = Ee^{e } W ] vf \ From (22) 
we obtain 

E Wa(q, a c )q 3 = ^(p + = (p + v)i , (23) 

3 i,j 

where we have made use of the primary constraints (15), in order to extend 
the sum on the RHS of (23) over all £, and of the completeness relation for 
the eigenvectors. Hence we have recovered (3) for a = a c . 

As we have pointed out above, the persistance of the primary constraints 
yields on Lagrange level equations involving only coordinates and velocities. 
These are part of the Euler-Lagrange equations of motion. From the point 
of view taken in this paper, that the equations of motion are obtained by a 
limiting procedure from those of an unconstrained system, the persistance of 
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the primary constraints can also be viewed to follow from the requirement 
that also the accelerations remain finite in the limit a — > a c . Thus for a ^ a c 
we have from (12) that <p£ = Xe q- v^ e °\ Taking the time derivative of this 
expression, and noting that lim Q ^ Qc \£ (q,a) = 0, and liuta,^^ di\e (q,a) = 
0, we immediately conclude that (pe (q,p,a c ) = 0. On the Hamiltonian level 
this requirement must be implemented explicitely to yield the missing Euler- 
Lagrange equations of motion, not manifest in (19) and (20). 

From the above discussion it is evident that the limit a — > a c must be 
carried out on the level of the Hamilton equations of motion of the uncon- 
strained system, whereas on the Lagrangean level we are allowed to take this 
limit directly in the Lagrangean. The equivalence between the Lagrangean 
and Hamiltonian formulations has been studied in detail in [3]. 

As an example consider the singular Lagrangean of the pure U(l) Maxwell 
theory: 

L\A\ >] = -- J d 3 xF^F^ . (24) 
Consider further the non-singular Lagrangean 

L[A», A»] = ~\j d 3 xF^F^ + -a J d 3 xA° , (25) 

which for a — * reduces to (24). This (non-covariant) choice is of course only 
the simplest one. Any other Lagrangean reducing to (24) in the appropriate 
limit would be just as acceptable. The canonical momenta conjugate to 
are given by 

n, = F^ + a5,o(d A ) . (26) 
The Lagrangean written in the analogous form to (1) is 

L = J d 3 xd 3 y \^(x,x°)W, u (x,y)A»(y,x°)\- J d 3 x V ^A(x))A^x)-V[A] , 

(27) 

where 

% := (0,-VA°) , (28) 
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and the matrix elements of the symmetrix matrix W read 

W 0l (x,y) = 0, 
W 00 (x,y) = aS {3 \x-y) , 
W^x,y) = 5^\x-y) . 

The potential V[A] is given by 



V[A]= fd*x (l Fij FV- l -{VA°Y 



A~ l] 2' 
and the canonical momenta, analogous to (3) read 

tt^x) = f d 3 yW, u (x,y)A»(y,x°) - 



(29) 
(30) 
(31) 



(32) 



(33) 



The matrix W possesses the following orthonormalized eigenvectors, labeled 
by a discrete and continuous index, replacing the discrete index I in (4) 



v^\x)=5 pu 5^\x-z) , 
and the corresponding eigenvalues of W are given by 

A<°>*> = a; = 1 (i = 1,2,3) . 

The Hamiltonian, analogous to (6) then takes the form 

H ^^|d 3 J|rf 3 x(7r M ( a :)+^(x))^(f) 



+ 



+ V[A] 



(34) 
(35) 

(36) 
(37) 



d 3 x(7T,(x)+^(x))^(f) 

Upon making use of (34), this expression reduces to 

H =\T.j d " z ^ + + \j + Vo) 2 + V[A] , (38) 

where in the present case r] = 0. The Hamilton equations of motion read 



A* = ^ = 7r i + dVL°, 

OTTi 

SH _ 
7Tj = -T- = . (39) 

For finite ^4° the first equation tells us that in the limit a — > 0, 7To must 
vanish, whereas ^4° remains completely arbitrary. Since 7To must vanish for 
arbitrary times, the third equation tells us that also V • 7? = 0. This is just 
the secondary (Gauss law) constraint. 

Alternatively we could have departed from a covariant form for an uncon- 
strained system by adding to the Lagrangean density in (24) the covariant 
term § {d^f. In this case r]^(A(x)) also depends on a, and following our 
general procedure one is led to the equations of motion 

A° = ^ = -n -V-A, (40) 

A* = ^ = 7r l + aM°, (41) 
dH 

in = -djF ji + <9V • (43) 

From the first equation it follows again that in the limit a — > 0, 7r = for all 
times, so that we are left with the standard equations for the pure Maxwell 
theory. We emphasize once more that the limit must be taken on the level of 
the equations of motion, and not in the Hamiltonian. The examples demon- 
strate the, of course, well known fact, that only after taking into account 
the primary constraints and their persistance in time, the full set of Euler- 
Lagrange equations of motion are generated. Thus primary constraints, and 
possible secondary constraints following from them, play a special role, while 
terciary, etc. constraints correspond to consistency relations hidden in the 
Euler-Lagrange equations of motion. With the primary constraints written 
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in the form (15), a strict iterative construction of the persistance equations 
for the constraints will necessarily parallel exactly the equations obtained on 
the Lagrangean level, irrespective of any possible bifurcations 4 . 

Let us summarize. In this paper we have shown that the Hamilton equa- 
tions of motion for any constrained system, in the form given by Dirac, can 
be obtained as a limit of the equations of motion for an unconstrained sys- 
tem. It was thereby shown that the primary constraints follow directly from 
the requirement that the velocities be finite in this limit. To obtain the 
full set of equations on Hamiltonian level, which translate into the Euler- 
Lagrange equations of motion, one must take into account the persistance 
of the primary constraints in time. These equations, which are implicit in 
the persistance requirement of the primary constraints, can also be viewed to 
follow from the requirement that also the accelerations remain finite in the 
above mentioned limit. The particular form of the term added to the singular 
Lagrangean which converts the system into a second class system, plays no 
role. The only requirement is that the unconstrained system reduces to the 
constrained theory of interest in an appropriate limit. We have demonstrated 
this for the case of the pure Maxwell theory. 
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